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摘要 :【 目 的] 研究 地 理 位 置 接近 的 山地 区 域 中 华 蜜蜂 Apis cerana cerana 的 形态 与 遗传 多 样 性 , 对 
开展 中 华 蜜蜂 遗传 资源 的 保护 具有 重要 意义 。 重 庆 多 山 , 如 大 巴山 武陵 山 .大 炎 山 等 ,当地 植被 茂 
盛 、 蜜 源 植物 丰富 ,适宜 中 华 蜜蜂 的 生存 繁衍 。 因 此 ,重庆 中 华 蜜蜂 是 研究 同一 地 域 复杂 地 形 下 中 
华 蜜蜂 形态 和 遗传 分 化 的 良好 对 象 。 本 研究 首次 全 面 调查 重庆 地 区 中 华 蜜蜂 种 质 资源 多 样 性 ,分 
析 不 同 生境 中 华 蜜 蜂 种 质 资源 的 分 化 ,考察 潜在 的 基因 交流 情况 ， SN 蜜蜂 的 资源 保护 提供 参 
考 。【 方 法 ] 运 用 吻 长 .起 长 .起 宽 、 起 面积 、 肘 脉 指数 、 第 3、4 背 板 长 等 形态 指标 对 来 自重 庆 不 同 山 
脉 和 周边 地 区 的 139 群 中 华 蜜蜂 进行 形态 测量 和 统计 ,同时 对 mtDNA ~ CO II E 3E 28 85 [X 
es tse 分 析 。【 结果】 形态 分 析 显 示 , 以 亚热带 季风 气候 为 主 的 武陵 山大 娄 山 
地 区 中 华 蜜蜂 能 分 以 亚热带 山地 气候 为 主 的 大 巴山 地 区 聚 群 ， 以 及 重庆 西北 部 的 聚 群 
Structure 群体 结构 uM 非 编 码 区 单 倍 型 序列 发 现 ， 不 由 四 二 二 六 中 华 蜜蜂 种 群 之 间 存 在 基 
因 交 流 ,而 两 大 山系 与 重庆 周边 地 区 中 华 蜜蜂 种 群 之 间 同 样 存在 基因 交流 。 此 外 发 现 了 3 个 新 的 
东方 蜜蜂 单 倍 型 ,分 别 是 大 巴山 CqH4、 大 类 山 CqH7 , AARAA CqH13。[【 结 de 
地 理 种 群 形态 发 生 了 分 化 ,表现 出 较 强 适应 性 。 各 山脉 和 地 区 中 华 蜜蜂 种 群 之 间 遗 传 结构 接近 , 具 
有 共同 遗传 单 倍 型 ,暗示 基因 交流 发 生 在 不 同 地 理 种 群 中 ,导致 重庆 中 华 蜜蜂 的 遗传 分 化 不 明显 。 
关键 词 : 中 华 蜜 蜂 ; 地 理 种 群 ; 形态 ; 遗传 多 样 性 ; 单 倍 型 ; 重庆 
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Morphological and genetic diversity of the Chinese honeybee, Apis cerana 














cerana ( Hymenoptera: Apidae) in Chongqing, southwestern China 

GU Ying' , BIN Xian-Li' , LIU Lu ,ZHOU Jun', MA Zhen-Gang' , TAN Hong-Wei^, ZHOU Ze-Yang! , 
XU Jin-Shan" ' (1. College of Life Sciences, Chongqing Normal University, Chongqing 401331 , China; 
2. Chongqing General Station of Animal Husbandry Technology Extension, Chongqing 401121, China) 

Abstract: [ Aim] Chongqing of southwestern China features hills and low mountains in which floral 
resource plants are abundant. As a consequence, there exists a rich resource of the Chinese honeybee, 
Apis cerana cerana in different geographical regions of Chongqing. Our study comprehensively surveyed 
the morphological and haplotype diversity of A. cerana cerana populations collected from different habitats 
in different regions of Chongqing, which will contribute to the resource conservation of A. cerana cerana. 
[Methods] Samples of one hundred and thirty-nine A. cerana cerana populations from different 
mountainous regions of Chongqing and surrounding areas were collected, and six morphological 
characteristics of these samples were measured. The fragments of mtDNA tRNA""-COJI non-coding region 
were amplified by PCR and then sequenced to evaluate their genetic diversity. [Results] Morphological 
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fi 瑛 等 : 重庆 中 华 蜜蜂 的 形态 和 遗传 多 样 性 








cluster analysis showed that honeybee populations from Wuling and Ta-lou mountains with subtropical 
monsoon climate could be distinguished from the honeybee populations from Ta-pa mountains with 
subtropical mountain climate, and the northern and western regions of Chongqing. The analysis of genetic 
structure of mtDNA non-coding region showed that gene flow occurred among the honeybee populations 
from different mountainous habitats. Three new haplotypes of A. cerana cerana were identified, i. e., 
CqH4 from Ta-pa Mountains, CqH7 from Ta-lou Mountains and CqH13 from western hills. [ Conclusion] 
The morphological characteristics of different geographical populations of A. cerana cerana in Chongqing 
show a degree of diversification and are consistent with the geographical distribution. 'The common 
haplotypes are usually shared between the A. cerana cerana populations from different mountains and 
regions of Chongqing, suggesting that gene flow occurs and thereby increases the genetic complexity of the 
Chinese honeybee. 


Key words: Apis cerana cerana; geographical population ; morphology; genetic diversity; haplotype; Chongqing 
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除了 提供 蜂 产 品 外 ,蜜蜂 在 维持 生态 系统 平衡 
和 促进 农作物 增产 方面 起 着 举足轻重 的 作用 ,例如 
蜜蜂 每 年 为 农作物 授粉 所 产生 的 经 济 价 值 高 达 数 百 
亿 元 (Sharmah et al.,2015 )。 中 华 蜜蜂 Apis cerana 
cerana 简称 "中 蜂 ”, 属 东方 蜜蜂 Apis cerana WPP, J 
泛 分 布 于 全 国 各 地 ,具有 悠久 的 饲养 历史 ( 陈 盛 禄 ， 
2001) 。 中 华 蜜蜂 可 利用 零星 蜜源 植物 , 采 蜜 期 长 ， 
非常 适合 中 国 山区 定点 饲养 ,但 自 20 世纪 初 西方 蜜 
蜂 引 入 以 来 ,由 于 西方 蜜蜂 产 蜜 量 大 ,中 华 蜜蜂 在 与 
西方 蜜蜂 的 种 间 竞 争 中 处 于 弱势 ,又 由 于 生态 环境 
遭 到 破坏 ,中 华 蜜蜂 资源 锐 减 ( 杨 冠 煌 , 2009) ,因此 
保护 中 华 蜜蜂 种 质 资 源 迫 在 眉 睫 。 

研究 中 华 蜜蜂 资源 分 化 是 进行 科学 保护 的 前 
提 , 目 前 对 蜜蜂 资源 分 化 研究 方法 主要 分 为 形态 学 
方法 和 分 子 生 物 学 方法 。 外 部 形态 是 物种 遗传 最 直 
接 的 表现 形式 ,利用 蜜蜂 的 形态 特征 标记 可 以 描述 
蜜蜂 多 样 性 (Ruttner et al., 1978; 李 志 勇 和 王 志 ， 
004)。 有 研究 者 较为 系统 地 对 西方 蜜蜂 的 形态 特 
征 进 行 了 鉴定 (Rutiner, 1988) , 目前 我 国学 者 对 于 
中 华 蜜蜂 多 采用 这 种 形态 分 析 方 法 。 人 研究 发 现 全 国 
各 地 区 的 中 华 蜜蜂 在 形态 上 发 生 了 分 化 ( 张 连 江 
等 , 2006; 王 桂 芝 等 ， 2008a，2008b; 罗 凌 绢 和 谭 
E, 2009; 朱 翔 杰 等 ,2009 ) ,海拔 、 经 纬度 、 年 均 气 
温 等 地 理 和 气候 因素 均 对 中 华 蜜蜂 分 化 有 着 一 定 的 
影响 ( 王 桂 芝 和 石 独 ，2008a，2008b; iE SE, 
2005 ) 。 在 分 子 水 平 上 对 蜜蜂 进行 遗传 分 化 分 析 是 
形态 学 鉴定 的 重要 补充 。mtDNA 标记 是 蜜蜂 遗传 
聚 类 的 主要 分 析 方 法 (Smith and Brown, 1988; 
Garnery et al., 1992, 1993; Smith et al., 2000; Frank 
et al., 2000) ,其 中 mtDNA tRNA™ ~ CO II 2A [8] B) 
一 段 富 含 AT 的 非 编 码 区 常用 来 进行 遗传 变异 研究 
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C amA R, 2005), Haiz E E 55 (2007 ) 利 
用 mtDNA 对 东方 蜜蜂 进行 了 分 析 , 发 现 11 个 地 理 
种 群 在 该 段 序 列 有 9 个 变异 位 点 ,揭示 了 东方 蜜蜂 不 
同 地 理 种 群 的 遗传 多 样 性 ;Magnus 等 (2014) 对 美国 
非 人 工 饲养 意大利 蜜蜂 种 群 mtDNA. 分 析 发 现 了 23 
个 单 倍 型 ,证 明 非 人 工 饲养 种 群 存在 着 不 同 的 谱系 来 
源 。 以 上 方法 为 我 们 研究 蜜蜂 资源 分 化 提供 了 参考 。 

重庆 位 于 长 江上 游 ,年 平均 气温 16 ~ 18Y , 平 
均 海 拔 400 m, 最 高 2797 m, 山 地 面积 占 70% ,有 大 
巴山 Ab RRL KALU. 复杂 的 山地 环境 
使 重庆 中 华 蜜蜂 资源 丰富 ( 曹 联 飞 等 , 2010; 于 增 
源 等 , 2012; 任 勤 等 , 2012; 谭 宏伟 等 , 2014 ) 。 
此 ,该 地 区 中 华 蜜蜂 是 研究 复杂 地 形 下 蜜蜂 形态 和 
遗传 分 化 的 良好 对 象 。 本 研究 通过 测定 形态 指标 测 
定 , 分 析 mtDNA 非 编 码 区 遗传 变异 与 群体 遗传 结 
构 , 对 重庆 不 同 山脉 及 不 同 地 点 的 中 华 蜜蜂 资源 多 
样 性 进行 了 调查 ,为 中 华 蜜蜂 资源 保护 提供 了 参考 。 
































1 材料 与 方法 


1.1 采样 

本 研究 共 采 集 139 群 中 华 蜜蜂 样本 ,主要 为 老 
桶 养殖 或 野生 招募 ,并 避免 在 有 引种 和 商业 育 王 的 
区 域 采样 。 样 本 采集 地 信息 见 表 1 ,其 中 重庆 共 采 
集 不 同 山脉 和 地 区 122 群 ,另外 采集 湖北 荆门 大 洪 
山 10 群 ,四 川 乐山 哦 眉山 6 群 ,阿坝 马尔 康 1 群 作 
为 比较 。 每 群 蜜蜂 采集 工蜂 50 头 左右 ,浸泡 于 
95% 乙醇 中 保存 。 
1.2 形态 指标 测定 

每 群 蜜蜂 选取 5 头 工蜂 , 按照 《中 国 畜 禽 遗 传 
资源 志 ( 蜜蜂 志 )》 所 要 求 测定 的 指标 内 容 及 以 往 相 
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R1 中 华 蜜蜂 不 同 地 理 种 群 采集 信息 
Table 1 Sampling data of different geographical populations of Apis cerana cerana 
"TT 采样 地 点 FERE 经 度 AE R) 
Mountains/ Regions Sampling locality Number of colonies sampled Longitude Latitude Altitude 
iX Wulong ( WL) 10 107°45'E 29°19'N 502 
武陵 山 Wuling Mountains 秀山 Xiushan (XS) 5 108°59'E 28?27'N 563 
HPH Youyang ( YC) 2 108?45'E 28*50'N 647 
I| Nanch NC 39 107*05'E 29?09'N 1 065 
大 娄 山 Ta-lou Mountains Mi) | E van CNE) 
基 江 Qijiang (QJ) 10 106?40'E 28*52'N 324 
庆 大 巴山 Ta-pa Mountains 城 口 Chengkou ( CK) 10 108?39'E 31?56'N 877 
* i KÆ Dazu (DZ) 7 105°44'E 29°39'N 407 
Chongqing 
合川 Hechuan ( HC) 5 106?16'E 29*58'N 242 
. 江津 Jiangjin (JJ) 9 106*15'E 29*17'N 850 
西北 部 地 区 D 9547 abor 
: 开县 Kaixian ( KX) 11 108?23'E 31?09'N 807 
Northern and western regions 
潼 南 Tongnan (TN) 1 105?48'E 30*08'N 436 
万 州 Wanzhou( WZ) 3 108?24'E 30?48'N 851 
永川 Yongchuan( YC) 10 105?55'E 29?2]'N 460 
湖北 eue 11" og2' 
Hubei 大 洪山 Dahong Mountains Jj] Jingmen (JM) 10 112?11'E 31?02'N 204 
ubei 
四 川 眉山 Mount Emei 乐山 Leshan (LS) 6 103?26'E 29?34'N 415 
Sichuan 阿坝 Aba( ABa) 1 101?42'E 32?54'N 2 600 


“括号 内 缩写 为 种 群 代码 Abbreviations in brackets are population codes. 





关 研 究 (国家 畜 禽 遗传 资源 委员 会 , 2011; RIE 
等 , 2012; 田 柳 青 等 , 2013 ) ,分 别 测定 每 头 蜜蜂 的 
吻 长 .右前 翅 长 . 宽 和 面积 \ 肘 脉 指数 .第 3 +4 HE 
长 6 个 指标 。 
1.3 mtDNA 遗传 多 样 性 分 析 

采用 动物 基因 组 DNA 快速 抽 提 试剂 盒 ( 上 海 
生 工 ) 提 取 中 华 蜜蜂 总 DNA ,保存 于 -20%C 备用。 
FH mtDNA tRNA ~ CO 工段 公开 引物 E2:5'-GGCA 
GAATAAGTGCATTG-3'; H2: 5'-CAATATCATTGATG 
ACC-3' ( Garnery et al., 1992) 进行 PCR f 3$, PCR 
体系 : 上 、 下 游 引物 各 1 pL (10 pmol/L), dNTP 
Mixture 4 uL (2. 5 mmol/L) , Taq fi 0. 25 uL (5 
unit/uL) , DNA 模板 2 pL (100 ng/uL) , ddH,O 
36.75 uL,10 x PCR buffer( Mg^* ) 5 uL, PCR 反应 
程序 :95% 5 min 预 变性 ; 94% 45 s 变性 ,50%C 45 s 
退火 ,72%C 1 min 延伸 , 共 35 个 循环 ;之 后 72% 延伸 
5 min , 胶 回收 后 进行 双向 测序 。 
1.4 数据 分 析 和 统计 

采用 SPSS 21 软件 对 形态 数据 进行 显著 性 分 
析 、 因 子 分 析 、 主 成 分 分 析 、 判 别 分 析 以 及 聚 类 分 析 
( Brannelly et al., 2015 ) 。 形态 指标 数据 采用 R 3. 13 
作 散 点 图 ( Chambers, 2008) , DNA 测序 结果 利用 
Mega 6 进行 多 重 序 列 对 比 (Tamura et al., 2013) , 
DnaSP 5 进行 单 倍 型 分 析 (Librado and Rozas, 
2009) , Arlequin 3. 5 进行 分 子 方差 分 析 ( Excoffier 
and Lischer, 2010) , Network 5 绘制 网 络 中 介 图 



































( Bandelt et al., 1999) , Æ NCBI 上 下 载 已 公布 的 东 
Jj 8E & mtDNA 非 编 码 区 单 倍 型 ,利用 DNAMAN 8 
进行 单 倍 型 序列 比 对 ,使 用 Structure 2.3.4 进行 群 
体 结 构 分 析 ( Pritehard et al., 2000) 。 





2 结果 


2.1 重庆 中 华 蜜蜂 不 同 种 群 的 形态 测定 
采集 的 139 群 中华 蜜 蜂 的 陶 长 (Pr) c BE. 
(FL) \ 右 前 翅 宽 (FB)、 右 前 怒 面 积 (AW)、 第 3 +4 
背 板 长 (T3 +4) 、 肘 脉 指 数 (Ci) 形态 指标 两 两 之 间 
的 线性 关系 如 图 1 所 示 。 通 过 显著 性 差异 水 平分 
析 , 表 明 来 自重 庆 武 陵 山 、 大 类 山 、 大 巴山 和 其 西北 
部 的 中 华 蜜蜂 在 吻 长 、 怒 长 、 翅 宽 、 翅 面积 .第 3 +4 
背 板 长 的 平均 值 均 大 于 来 自 大 洪山 (湖北 研 门 ) 和 
峨眉 山 ( 四 川 乐山 ) 的 中 华 蜜蜂 。 其 中 大 巴山 (重庆 
城 口 ) 中 华 蜜 蜂 的 形态 指标 优势 明显 ,接近 于 采集 
自 阿坝 马尔 康 的 中 华 蜜蜂 。 
2.2 重庆 中 华 蜜蜂 种 群 间 形 态 指 标 判别 分 析 
利用 逐步 判别 分 析 的 Wilk”? Lambd 值 对 6 项 指 
标 进 行 分 类 ,绘制 判别 函数 散 点 图 。 根 据 判 别 质心 
将 不 同 地 理 中 华 密 蜂 种 群 进行 聚 类 ,判断 不 同 地 理 
种 群 中 华 蜜蜂 的 形态 距离 。 结 果 如 图 2 所 示 ,来自 
于 重庆 武陵 山 和 大 类 山 以 及 西 .北部 的 中 华 蜜蜂 至 
为 一 个 整体 群 ,而 来 自 于 重庆 大 巴山 的 中 华 蜜蜂 与 
上 述 整 体 群 较 远 。 阿 坝 马 尔 康 种 群 与 其 他 中 华 蜜蜂 
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图 1 重庆 中 华 蜜蜂 不 同 种 群 间 6 项 形态 指标 矩阵 散 点 分 布 情况 


Fig. 1 The matrix scatters of six morphological indexes of Apis cerana cerana populations in Chongqing 
Pr; 吻 长 Length of proboscis; FL: 右前 翅 长 Length of right forewing; FB: 右前 翅 宽 Width of right forewing; AW: 右前 翅 面 积 Area of right 
forewing; T3 +4: 第 3 +4 1$ c I& Length of tergum 3 +4; Ci; 肘 脉 指数 Cubital index. 





Canonical discriminant function 





1 武陵 山 Wuling Mountains 

2 大 娄 山 Ta-lou Mountains 

3 大 巴山 Ta-pa Mountains 

4 西 、 北 部 地 区 Northern and western. regions 
5 大 洪山 Dahong Mountains 

6 阿坝 ABa 

7 峨眉 山 Mount Emei 

组 质心 Group centroid 
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Function 2 








Function 1 


图 2 重庆 中 华 蜜蜂 不 同 种 群 间 6 个 形态 指标 判别 分 析 
Fig. 2 Discriminant analysis of six morphological indexes of Apis cerana cerana populations in Chongqing 
每 个 圆圈 代表 来 自 某 山脉 或 地 区 的 一 个 种 群 ,圆圈 颜色 代表 不 同 山 脉 或 地 区 ; 蓝 色 方 框 代表 每 组 质心 。 Colorful circles represent the honeybee 


populations derived from different geographical regions. Blue square represents centroid of each group. 
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表现 出 明显 差异 ,但 与 大 巴山 距离 相对 较 近 。 大 洪 

L CHAE) 与 峨眉 山 ( 四 川 乐山 ) 中 华 密 蜂 距 离 

较 近 ,而 与 其 他 地 区 较 远 。 

2.3 重庆 中 华 蜜蜂 种 群 间 形态 指标 聚 类 分 析 
进一步 利用 平方 欧式 距离 法 对 形态 指标 进行 直 

接 聚 类 (图 3) 。 结 果 显 示 , 来 自 于 重庆 不 同 山 脉 和 

地 区 的 中 华 蜜蜂 与 大 洪山 、 峨 眉山 的 中 华 蜜蜂 明显 























分 化 为 两 支 ,与 判别 分 析 结 果 相 一 致 ;而 重庆 中 蜂 主 
要 分 为 三 大 文 , 第 1 LEAR KAUR, ERA 
武 隆 、 南 川 等 ;第 2 文 是 大 巴山 聚 群 ,主要 以 城 口 为 代 
表 ; 第 3 支 是 以 合川 为 代表 的 西 .北部 聚 群 。 这 表明 
中 华 蜜蜂 在 不 同 地 理 条 件 下 形态 出 现 了 分 化 。 来 自 
于 武陵 山 和 大 娄 山 的 中 华 蜜蜂 聚 群 较 近 ,这 与 武陵 山 
和 大 类 山 均 属于 亚热带 季风 气候 的 地 理 环境 相关 。 
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图 3 基于 6 个 形态 指标 的 

















P 华 蜜蜂 不 同 地 理 种 群 形态 聚 类 








Fig. 3 The clusters of different geographical populations of Apis cerana cerana based on six morphological indexes 
种 群 代码 同 表 1, Population codes same as Table 1. 


2.4 mtDNA 非 编 码 区 的 碱 基 组 成 及 变异 情况 

测序 获得 了 139 群 中 华 蜜蜂 mtDNA tRNA™ ~ 
CO 工段 的 部 分 序列 ， 其 中 非 编 码 区 序列 大 小 为 97 
bp,A +T 含 量 为 87.29% ,是 一 段 A +T 富 集 区 。 对 
mtDNA JETI DX EZ B2 Ze FE RE ETT 23 PT ,结果 显示 
重庆 中 华 蜜蜂 的 平均 核 昔 酸 差异 数 (K) 29 1.7021, 
TECH EZ RERECPI) X 0.01792, FER EFIA 
RERA E RE RERED REMA E E E o 

分 析 中 华 蜜蜂 mtDNA 非 编码 区 序列 , 共 发 现 了 
17 个 单 倍 型 。 与 从 NCBI E P 2X BS ZR 7; 3E E 
mtDNA 非 编 码 区 单 倍 型 进行 对 比 ,发 现 了 3 个 新 的 
单 倍 型 (图 4) ,分 别 命名 为 CqH4 型 (大 巴山 )、 
CqH7 型 (大 类 山 ) ,CqH13 型 (西部 丘陵 ) 。 
2.5 mtDNA 非 编码 区 的 单 倍 型 网 络 中 介 分 析 

将 获得 的 mtDNA 非 编码 区 单 倍 型 变异 位 点 绘 
























































制 网 络 中介 图 ,结果 显示 CqH1, CqH3 和 CqH5 为 
3 个 主体 单 倍 型 单 倍 型 ,CqHl 和 CqH5 出 现 的 频率 
最 高 ,包含 的 地 区 最 丰富 ,为 主要 的 共有 单 倍 型 (图 
5)。 其 中 有 10 个 单 倍 型 起 源 于 CqH1,1 个 单 倍 型 
起 源 于 CqH5, CqH7 与 其 他 单 倍 型 表现 出 分 离 。 比 
较 来 自重 庆 不 同 山脉 和 地 区 的 与 大 洪山 阿坝 马尔 
康 以 及 峨眉 山 的 中 蜂 单 倍 型 时 发 现 ,它们 之 间 存 在 
共有 单 倍 型 ,并 没有 出 现 分 离 。 
2.6 mtDNA 非 编 码 区 群体 结构 分 析 

基于 形态 指标 聚 类 结果 将 样本 划分 为 5 个 地 理 
育 群 ,对 测序 获得 的 线粒体 序列 进行 群体 结构 分 析 。 
结果 如 图 6 所 示 ,将 等 位 基因 组 分 数 和 矩阵 种 群 数量 
(K) 设 置 为 2, 3 和 4 时 ,所 调查 的 5 p rp AERE MERE 
群 的 遗传 结构 整体 差异 性 均 不 明显 ,表明 来 自重 庆 
不 同 山脉 和 地 区 的 中 蜂 种 群 之 间 以 及 与 大 洪山 和 路 
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Hunzhou TGT GTC TGC GA 

CqH1 === == = UH 6 na 
CqH2 -一 一 下 -A- —- —T —  Yunnan2 一 -一 -一 
Cais, == = T= CH 8 be ne 
CqH4 -A- -— — G H9 = lil Se 二 
Cglb | == = TAT-- CH 10 coc /== Gansu —— BERG 
CqH6 -—- —- kc- — C€OH1l cn Shanxi2 | -—— ———— — 
CqH7 --C —TA4A---- CH 12 A nly 
CqH8  --- --- A-T -- CH 13 aA = = Hanan M ESSI 
CqH9 = = ,Jiling DUMP em eee 
CaH10 -— -— Me -- CH 15 eu Se -- dquüangdong ——— -—- -9—- se 
CqHll 一 一 A — — CHI 16 pA-es eec dambodialo-—--—Ac- i A- 
Coal2 è —— A— —- — CH 17  — —- —- —- —  Cambodia2 -A- 一 一 一 
CqH13 人 er Vietnam === (o —— 
CqH14 -—- -——— —-T —— CH 19 => Ana Veta e 生生 
CqH15 -— — -A- A- Gansu nl 
CqH16 -A- —- ——- — Shanxi r ==  Daosl cdm cM 
CqH17 ==> === == dojehüano  —— -—— cA- —— Maisl = TESA 
Cqne Hl -一 -一 -一 一 Henan 二 GIC 
Cqnc H2 —- -—— —- --  Nepall —— e oes Thail? rax a LE 
Cune Hà. == 4 ——- A-  Nepal2 C—eoeesscc Indias Hes 59 A -—- -—— 
Cqne l4 — —— A —- A- Nepal | —- A— -A- — India H3 一 一 T AT 一 
CH 1 -—- -—— —T-- Nepal4  —- —- —- — IndiaH2 一 一 A —- A- 
CH 2 -一 一 -一 Qinghai —- A— —- — Jaan 一 一 一 一 
CH 3 一 -一 -并 一 Sichuani 一 -一 一 一 mdiap2 -A- A-T -TA TT 
CH 4 -—- -—— —T -—- Sichuan2 一 - -一 -一 一 IndiaPl -A- A-T -TA TT 
CH 5 Teca ehian v 9e 

到 4 ”东方 蜜蜂 不 同 种 群 间 mtDNA 非 编码 区 单 倍 型 序列 对 比分 析 














Fig. 4 Alignment of haplotypes of mtDNA non-coding region in Apis cerana populations 
CqHI1 - CqH17 代表 本 研究 所 测定 获得 的 17 个 单 倍 型 ,红色 标记 表示 新 发 现 的 单 倍 型 。 短 横 线 代表 相同 序列 。CqHI - CqH17 are the symbols 


of seventeen haplotypes sequenced in this study. Symbols marked with red represent the haplotypes not ever identified. The short lines mean the common 


nucleic acids. 
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图 5 鉴定 获得 的 重庆 中 华 蜜蜂 不 同 种 群 mtDNA 非 编码 区 单 倍 型 网 络 中 介 分 析 


Fig. 5 The median-joining networks with the haplotypes of mtDNA non-coding region from Apis cerana cerana populations in Chongqing 





区 | 








中 每 个 圆圈 代表 一 种 单 倍 型 ,其 面积 大 小 与 该 种 和 














倍 型 出 现 的 频率 成 了 





E 比 。 红 色 数 字 代表 两 两 





feo pa] AE Su. Each circle represents 





one haplotype and its area stands for the frequency of haplotype occurrence. The numbers marked with red denote the site of mutation between two 


haplotypes. 
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图 6 重庆 不 同 地 理 聚 群 的 中 华 蜜蜂 mDNA 非 编 码 区 单 倍 型 的 遗传 结构 组 成 


Fig. 6 The genetic structure of haplotypes of mtDNA non-coding region based on morphological cluster of 


Apis cerana cerana populations in Chongqing 





K 代表 种 群 数量 。K represents the number of populations. 


眉山 的 中 华 蜜蜂 种 群 虽然 在 形态 上 有 差别 ,但 遗传 结 
构 组 成 相似 。 这 一 结果 暗示 了 基因 交流 广泛 发 生 在 
不 同 地 理 分 布 的 中 华 蜜蜂 之 间 ,使 得 遗传 分 化 与 地 理 
分 布 差异 性 乃至 形态 多 样 性 之 间 的 关联 性 非常 小 。 











3 讨论 


蜜蜂 是 重要 的 经 济 性 授粉 昆虫 ,研究 中 华 蜜蜂 
形态 和 遗传 分 化 可 以 了 解 我 国 不 同 地 区 的 蜂 种 资 
源 , 为 中 蜂 的 生态 资源 保护 提供 指导 。 本 人 研究 采集 
了 重庆 不 同 生 境 122 群 中 华 蜜蜂 样本 ,测量 和 分 析 
Y 6 项 形态 指标 ,并 与 大 洪山 (湖北 荆门 ) ,马尔 康 
(四 川 阿 坝 ) 、 峨 眉山 (四 川 乐山 ) 的 17 群 中 华 蜜蜂 
形态 指标 进行 比较 。 结 果 显 示 重 庆 不 同 山脉 和 地 区 
的 中 华 蜜蜂 吻 长 、 怒 长 翅 宽 以 及 第 3 和 4 背 板 之 和 
整体 有 别 于 大 洪山 、 峨 眉山 中 华 蜜蜂 。 由 于 重庆 地 
形 复杂 ,海拔 落差 较 大 ,形态 上 的 差异 应 该 与 生境 适 
应 相关 ( 曾 志 将 ,1992; Hepburn et al., 2000; RAK 
生 等 , 2010) 。 直 接 聚 类 和 主 成 分 聚 类 均 显 示 重 庆 中 
蜂 与 大 洪山 中 蜂 明 显 分 化 为 两 大 类 , 这 暗示 着 重庆 的 
中 华 蜜蜂 可 能 与 华中 型 中 蜂 在 形态 上 有 区 别 。 本 研 
究 所 获得 的 重庆 中 蜂 3 个 形态 聚 群 中 ,大 巴山 为 北 亚 
热带 山地 气候 ,其 他 地 区 均 为 亚热带 季风 气候 ,气候 
类 型 和 所 处 山脉 的 不 同 使 中 华 蜜蜂 形态 产生 分 化 ,这 
































与 前 人 的 研究 结果 相 一 致 ( 余 林 生 等 , 2008) 。 

蜜蜂 具有 独特 的 婚 飞 繁殖 方式 ,工蜂 的 mtDNA 
能 反映 蜂 群 的 遗传 结构 ,因此 mtDNA 标记 是 研究 蜜 
蜂 选 传 多 样 性 的 依据 ( Crozier and Crozier, 1993 ) 。 
不 同 地 理 种 群 中 蜂 mtDNA 存在 差异 ( 刘 振 国 等 ， 
2015) ;本 研究 发 现 ,同一 地 点 的 不 同样 本 其 mtDNA 
序列 也 存在 差异 ,如 大 巴山 脉 的 城 口 就 有 4 个 单 倍 
型 。 本 研究 鉴定 了 17 个 单 倍 型 ,对 比 发 现 有 3 个 新 
的 单 倍 型 ,分 别 来 自 大 巴山 脉 ( 城 口 ) 、 大 类 山脉 ( 南 
川 金 佛山 ) 、 西 部 丘陵 (合川 ) ,这 表明 中 华 蜜蜂 在 重 
庆 复 杂 地 形 下 种 群 遗传 资源 丰富 ,具有 很 大 的 保 种 
价值 。 重庆 不 同 山 脉 和 地 区 的 中 华 蜜 蜂 之 间 
mtDNA 遗传 结构 类 似 ,与 大 洪山 和 峨眉 山 的 中 华 蜜 
蜂 也 相 类 似 。 可 见 , 地 理 环境 造成 的 生境 条 件 差 异 
并 不 是 影响 中 华 蜜蜂 种 群 间 分 子 亲 缘 关 系 的 决定 因 
素 。 不 同 生境 条 件 下 蜜蜂 种 群 之 间 的 基因 交流 会 影 
响 中 华 蜜蜂 的 遗传 分 化 ( Bertrand et al., 2015, HE 
萍 等 , 2016) ,这 是 导致 重庆 中 蜂 不 同 地 理 种 群 之 间 
遗传 结构 差异 性 不 明显 的 原因 之 一 。 此 外 ,形态 分 
化 受 多 因素 影响 , 如 核 基因 的 转录 水 平 差异 .DNA 
表 观 遗传 调控 等 ,而 mtDNA 是 母系 遗传 分 子 ,主要 
反映 系统 遗传 分 化 。 因 此 ,重庆 中 华 蜜蜂 不 同 地 理 
种 群 的 形态 变化 与 mtDNA 遗传 分 化 的 关联 性 不 
明显 。 
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